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Abstract 

The fractal dimensions of the hull, the external perimeter and of the red bonds are measured through 
Monte Carlo simulations for dilute minimal models, and compared with predictions from conformal 
field theory and SLE methods. The dilute models used are those first introduced by Nienhuis. Their 
loop fugacity is |3 = — 2cos{7t/K) where the parameter i< is linked to their description through confor- 
mal loop ensembles. It is also linked to conformal field theories through their central charges c{R) = 
13 — 6(k + K^') and, for the minimal models of interest here, k — p/p' where p and p' are two co- 
prime integers. The geometric exponents of the hull and external perimeter are studied for the pairs 
(p,p') = (1,1),(2,3),(3,4),(4,5),(5,6),(5,7), and that of the red bonds for (p,p') = (3,4). Monte Carlo 
upgrades are proposed for these models as well as several techniques to improve their speeds. The mea- 
sured fractal dimensions are obtained by extrapolation on the lattice size H, V — > oo. The extrapolating 
curves have large slopes; despite these, the measured dimensions coincide with theoretical predictions up 
to three or four digits. In some cases, the theoretical values lie slightly outside the confidence intervals; 
explanations of these small discrepancies are proposed. 

Keywords: dilute logarithmic minimal models; logarithmic minimal models; conformal field theory; con- 
formal loop ensembles; SLE; fractal dimensions; geometric exponents; Monte Carlo simulations. 
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1 Introduction 



Geometric objects remain a central tool in the study of the critical behavior of statistical lattice models. 
Some of the most natural ones are the hull of a spin cluster, its mass, external perimeter and red bonds. 
Even though they were investigated as early as in the 1970's, their role has remained central over the 
years. In the 1980's, their close ties to conformal weights of the Virasoro algebra, and more generally 
with conformal field theory (CFT), was revealed starting with works by Saleur and Duplantier 1241 . And 
in the late 1990's, techniques from probability theory related their properties to that of random curves 
grown through stochastic Loewner evolution (SLE). The goal of this paper is to measure, using Monte 
Carlo simulations, the fractal dimensions of these objects for dilute lattice models. 

Recently, Sarnt-Aubrn, Pearce, and Rasmussen Il23l measured these dimensions for a family of lattice 
loop models whose continuum scaling limit is called the logarithmic minimal models (Pearce, Rasmussen, and 
Zuber ||20] ). Their results gave compelling evidence for the theoretical predictions of Saleur and Duplantier 
Il24l and others, and confirmed the rigorous result by Beffara IIH for the hull fractal dimension. Their work 
probed the dense phase of loop models and the present paper can be seen as completing their work by 
considering the dilute phase. 

We shall do so on another family of loop models based upon the celebrated ©(n) model. Writing the 
loop fugacity as n = —2 cos these loop models are well defined for all real k values, and our methods 
apply for arbitrary values of k. Focus here is on rational values, though, for which the fractal dimensions 
are rational and expressible in terms of conformal dimensions from an extended Kac table. We furthermore 
believe that, for k rational, these loop models converge to logarithmic CFTs in the continuum scaling limit 
and henceforth refer to them as dilute logarithmic minimal models. 

Beside the intrinsic value of checking theoretical predictions through experiments or, in the present 
case, Monte Carlo simulations, such checks often lead to improvements in the techniques of (numerical) 
experimentation. Together with Potts models, the XXZ Hamiltonian and other spin lattice models, dilute 
loop models are of great theoretical value. But, because the Boltzmann weights of loop models require the 
coiinting of the number of loops in a configuration, a task that is highly non-local, simulations of these loop 
models remain difficult. The classical algorithms, like that of Swendsen and Wang ||27|| , usually do not apply 
to them. There has been progress to extend these cluster algorithms to larger families of models, e.g. by 
Chayes and Machta IZl and Deng, Garoni, Guo, Blote, and Sokal ||8l but, unfortunately, some loop models 
remain without efficient algorithms. One of the outcomes of the present paper is the proposed upgrade 
algorithm and its variants that curtail significantly the difficulties of visiting large non-local objects, like 
loops whose size is commensurate to that of the lattice. 

This paper is organized as follows. The next section recalls the definition of the O (n) model and char- 
acterizes the interval of the parameter that will lead to their dilute phase. The definition of the geometric 
fractal dimensions for the hull, external perimeter and red bonds are given with their theoretical predic- 
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tions. Section 3 describes difficulties inherent to measuring geometric exponents on dilute models: the area 
associated to the various states of boxes, the important boundary effects and the distribution of the defect's 
winding number. Section 4 presents the results and discusses some of their shortcomings, while section 5 
contains some concluding remarks. Technical details are gathered in the appendices: the Monte Carlo al- 
gorithms are described in appendix A, and the statistical analysis in appendix B. Appendix C recalls the 
definition of the dense loop models used in 12311 . 



2 The 0[n] model and dilute logarithmic minimal models 
2.1 Loop representation of the two-dimensional 0{n) model 

The two-dimensional 0[n) model is a lattice spin model, whose n-component spins s^, for |x e {1 , 2, • • • , n}, 
are located at each site i G 5 of the lattice, with S being the set of sites. The spins are constrained to live on 
an (n — 1 )-dimensional sphere of radius y/n, i.e. they satisfy s, • Si = ^^=i s^s^ = ^i^t ~ ^- indicated, 
no summation is implied on Latin letters. 

2.1.1 On the honeycomb lattice 

The partition function of this model on a domain of the honeycomb lattice is defined as Z = Tr 0(ij> ^"''^ 
where the product is taken over all lattice edges, written (ij), and x is the inverse temperature. (See for ex- 
ample Dubail, Jacobsen, and Saleur IfTOl .) The high-temperature expansion, that is for x small, is given 
by 



Tr 



The trace operator 



Tr[A(si„Si„...)] =cir[f 6(s['sr-T^)dSiiA(Si„Si„, 



(1) 



(2) 



with C a normalization factor, allows to compute the expectation value of an arbitrary function A of the 
spins. It may be normalized such that the following properties hold: 



TrH] = 1 

Tr[srs^]=6,,6^- 
Tr[sr]=Tr[(sO' 



Tr (s^y 



nn 

2r(^ + V 



(3) 



with 151 the number of sites in the 



In this case, the normalization factor is given by C = |^ 
domain. With these properties, the partition fimction ([T]| becomes a sum on all configurations with non- 
intersecting loops. This is so because only cyclic terms of the form x^s^' s[^' s^^s 



g^i^k |ik sui-vive the 

l-k 1-1 
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trace, yielding a weight nx^ for the configuration. Applying these rules leads to the celebrated loop parti- 
tion function 

c 

for the 0[n) model on the honeycomb lattice. The sum is taken over all lattice configurations £ of non- 
rntersecting loops. Here X and N are respectively the total number of monomers (bonds) and the total 
number of loops of the configuration. The parameters x and n play accordingly the roles of bond and loop 
fugacity. The loop representation (|4| describes a larger family of models than the spin representation 
because n may take real values, not only positive integer ones. 

2.1.2 On the square lattice 

The trace operator ||2} forbids the presence of terms of odd powers in s^, so diagrams of intersecting or 
open loops are impossible in any C on the honeycomb lattice. But loop intersections may occur in the 0[n) 
model ([ij on the square lattice, because terms like x^[s^)^ are possible there and non-zero in general. The 
square-lattice O (n) model we are interested in is not equivalent to although it is quite similar, and was 
first defined by Nienhuis Ill9l . This model is defined on the square lattice, possesses a partition function 
similar to Q and has no intersecting loops. 

The spins of the model are located on the edges of the square lattice and there is now three interaction 
constants: u, v, and w. These parameters are understood to include the inverse temperature. As for the 
previous model, the spins are constrained by s[^s[^ —n and the partition fimction is defined as 



Z = Tr 



Y[ Q(Si.Sj,Sk,Sl] 



(5) 



where the product is taken over all lattice faces (i, j, k, I) with i, j, k and I the spins surrounding the face in 
fixed order (say clockwise). The Boltzmann weight Q is of the form 

Q(Si,Sj,Sk,Si) = 1 +u(Si • Sj +Sj • Sk + Sk -Si + Si • Si) +v(Si • Sk + Sj - Si) +w[(Si -SjJfSk • Si) + (Sj • Sk)(Si - Si)]. 

Note that all possible pairings of the face spins appear except (s, • Sk)(Sj • Si), which would allow for inter- 
secting loops. The trace operator Tr in ^ is defined again by Note that the model is given here at the 
isotropic point, a simplification first given by Blote and Nienhuis IS). 

Using the properties l|3| of the trace, the partition function ||5} becomes 

Z^^u^-v^-w^-n'^. (6) 

As before, £ is the set of non-intersecting loop configurations and n^, and n^, are respectively the 
numbers of u, v and w type faces (see figure [ij in the configuration, and N is still the total number of 
loops in C. Again, n is the loop fugacity while u, v and w are the weights of each face type. In Blote and 
Nienhuis fS], the set of critical points is divided into five non-equivalent branches parametrized by n. We 
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Figure 1: The Boltzmann weight of the nine possible faces. 



are interested in two of them here, corresponding respectively to the dilute and dense phases. They are 
both parametrized by 



sin 2A sin ^ 

u — 



sin 2A sin 3A + sin^ ^ 



2 3A 

(7) 



W : 



sin 2A sin 3A + sin^ ^ 

sin I sin ^ 
sin 2A sin 3A + sin^ ^ 



n = — 2 cos 4A, 

with the crossing parameter A constrained by < A < J or ^ < A < J for the dilute and dense phases, 
respectively. Note that, in both phases, the parameter n covers the range [—2, 2] once. The weight of the 
empty face is normalized to 1 for all A. The original weights given in Blote and Nienhuis ||5| can be obtained 
by replacing A with ^ (tt — 9) in ||7|. 

We conclude this subsection by recalling the full weights of the dilute loop models ||T91 by including the 
spectral parameter here indicated by ft to distinguish it from the weight u that appears in Q and figure [T| 
The weight of the empty face is 

sin-Clsin(3A — ■Q.l 

1 H ^ 

sin 2A sin 3A 



while 



sinU sin(3A — H) 

sm 3A sm 3A 

sin-Q.sin(3A — H) 

^ = ■ -,1. ■ 11. — (8) 

sm 2A sm 3A 

sin{lsin(A--a) sin(2A - -a) sin(3A - 11) 

W] = , W2 = 

sin 2A sin 3A sin 2A sin 3A 
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where the weights of the faces in the second and fourth columns in figure [T| are labeled from below by 
U] , U2, U3, U4 and W] , W2, respectively. At the isotropic point H — 3A/2, these weights reduce to the ones in 
||7| after rescaling them to get the weight 1 for the empty face. 

2.2 Dilute logarithmic minimal models 

The descriptions of physical systems through spin. Potts and loop models usually have different transfer 
matrices. One striking difference is the dimensions of the vector spaces upon which the spin or loop trans- 
fer matrices act; these dimensions are not equal. It is therefore not surprising that different continuum 
scaling limits may coexist for the same model, depending on the spin or loop description under study. By 
studying a dense loop representation of a family of models, different from the one used here, Pearce et al. 
II20I found that some of the associated transfer matrices exhibit nontrivial Jordan blocks. They subsequently 
argued that this gives rise to logarithmic CFTs in the continuum scaling limit and labeled a two-parameter 
family of such limits by CA4(p,p'), for logarithmic minimal models, where p and p' are positive coprime 
integers (see appendix[C||. Similarly, we believe that logaritmic CFTs arise in the continuum scaling limit of 
the loop models defined by ||6]|-||8]| and thus label the corresponding two-parameter family of continuum 
scaling limits as 2?£A^(p,p'), for dilute logarithmic minimal models, where p and p' are as above. One way 
to justify the prefix "dilute" is by the visual aspects of the loop configurations, as opposed to those of log- 
arithmic minimal models in which only w-type faces are admissible. See figure |2]f or typical configurations 
of 2?£A^(p,p') and appendix |C] for one of £A^(p,p'). 

Because loops in both the dense and dilute phases of X'£7W(p,p') are likely to be related to conformal 
loop ensembles CLEk (see Camia and Newman Werner ESlI , Sheffield 1251 ), the relationship between the 
pair (p, p '), the parameter k (or k = k/4) and the crossing parameter A needs to be given. The (logarithmic) 
CFT underlying I?£Al(p,p') has central charge 

/ I \ I ^ in the dilute phase 

c = c(k) = 13 - 6 K + - , whereR=< , (9) 



K 



in the dense phase 



and conformal weights 



A,,, =A,,s(k) = ^ '—^ r,s = 1,2,3,... (10) 



4k 

The usual duality k ^ of CLEk here becomes i< ^ and implies 

0=c(k) and A,,s =As,r(K). (11) 
The link between 2?£A^(p,p') and the 0[n) model is completed by the expression of A in terms of k: 

A(.) = ^ (l - ^) • (12) 

This form for A( R) is due to our choice of parametrization in |[7|; one can also find A( R] — ^ in the literature. 
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In terms of k, the dilute and dense branches correspond to 

^ < K < 1 , dilute phase, 

2 (13) 
1 < R < oo, dense phase. 

Moreover, we shall henceforth rename the loop gas fugacity n by (3: 

(3 =-2cos4A = -2cos(^) . (14) 

The duality R i and its implication < (TT) suggest that the dense and dilute phases of I?£A^(p,p') 
are dual to each other. There also exists a link between 2?£Al(p,p') and £A1(p,p'), which is twofold. To 
appreciate these dualities and links, we say that two models belong to the same universality class if they 
have the same central charge and share the same Kac table of conformal weights Ar,s (i<); this is typical of 
what is found in the CFT literature. (At this point, we do not require, or even address, that the models 
are based on the same set of representations with identical Jordan-block structures.) First, it is easy to see 
that the dense phase of VCAi (p , p ' ) lies in the same universality class as CAi (p , p ' ) since, for a given pair 
(p, p '), both models have the same R (eq. l|9|), and thus the same central charge c and Kac table. Second, the 
dilute phase of VCM. (p , p ' ) belongs to the same universality class as CM (p ' , p ) . This is also obvious since, 
by inverting p p ' in the definition of R in the logarithmic minimal models, one falls back to the R of the 
dilute phase of VCM. (p , p ' ] . Consequently, the conformal weights and the central charges are also equal in 
these models. 

This equivalence has been known for a while: Duplantier IITTl [121 used the Coulomb gas picture to 
discover a correspondence between the loop representation of the Ofn) model and the Fortuin-Kasteleyn 
(FK) representation of the -Potts model, which is valid for both the dilute and dense phases, when n is 
in the physical regime, that is n e [0, 2]. Since CM is built upon the FK representation of the Potts model 
and VCM on the O (n) model, the equivalence merely appears by inheritance. 

The principal series, i.e. models with p ' = p + 1 , of the dilute logarithmic minimal models are 

dense phase of 2?£A^(p,p') : 

i<dense = ^ = {2) f , |) |> f i • • • , 1 } <^ {polymers, percolation, Ising, tricritical Ising, 3-Potts, . . . , 4-Potts} 

(15) 

dilute phase ofI?£A^(p,p') : 

Kdiiute — fy — {i)f>f>t)|>---)l} ^ {polymers, percolation, Ising, tricritical Ising, 3-Potts, . . . , 4-Potts}. 

In the case of VCMll, 1), p and p' are strictly speaking not coprime. However, as indicated in | (15) , the 
model can be viewed as arising from the limit p ^ oo of the principal series, for which p and p' are 
coprime. For our purposes, though, it suffices to simply set p = p ' = 1 . 

The models in (TS) are named according to their universality class. For instance, it is well known that 
the Ising model belongs to the universality class of c = 1/2; this is why the models with Rdense — 4/3 
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(Pdense — V2) and Kdiiute = 3/4 (pdiiute — 1) are both called "Ising model" here. This might be confusing 
however since, in the 0(rL] model literature, one reserves that name for 0{n — ^ — ]), for both the dilute 
and dense phases. For comparison, a typical configuration for both phases of 2?£A^(3,4) is shown in 
figure |2] 





(a) Dense phase configuration. 



(b) Dilute phase configuration. 



Figure 2: Typical configurations on the 32 x 32 cylinder for the dense and dilute phases of VCM (3, 4) with 
Pdense — V^, pdiiute = 1/ and c = 1 /2 for both. Note the large loop in the dilute configuration. In both 
phases, such large loops often occur. 



2.3 Geometric fractal dimensions 

The hull perimeter, or simply hull, is defined as the set of all outer boundary sites of a cluster (see Grossman 
and Aharony IITSl ). Since here there are no clusters, but only loops, this interpretation has to be adapted to 
the present models. To do this, we use the same trick as in 1231 where it was applied to the CAi models. We 
thus consider that the loops present in VCAd configurations are Peierls contours of clusters of spins living 
on the lattice sites. Moreover we shall choose the boundary conditions such that at least one loop enters 
one of the boxes in the top row and exits in one of the bottom rows. This "loop" must then cross vertically 
the whole lattice and measurements can be done on this single object instead of taking averages over loops 
in each configuration. Such a loop will be called a defect. The analogy with the spin definition can then be 
constructed easily. The hull for a configuration with two defects is represented in pale (yellow) in figure |3] 
The external perimeter is the set of accessible outer boundary sites of a cluster and constitutes a subset of 
the hull. In loop models, the external perimeter of a given loop can be interpreted as its contour minus 
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Figure 3: The three observables. This figure was created specifically for the purpose of illustration and not 
randomly generated for a particular model. 

the set of fjords. For dilute models, a fjord is created each time the defect bounces on itself by using a w- 
type face, entrapping an area which now becomes inaccessible to the defect. (We have borrowed the term 
fjord from Asikainen, Aharony Mandelbrot, Rauch, and Hovi fS] where fjords with various gate sizes are 
considered. Ours have the narrowest gate size possible.) This is the discrete analogue of the situation where 
an SLE path engulfs a domain in its hull. It is emphasized that the definition of "hull" used here is different 
from the one used in SLE. The external perimeter is shown by a darker (green) line in figure |3| sitting above 
the paler hull. 

The red bonds, also known as singly connected bonds, were first introduced by Stanley 126] . They are 
defined as the bonds of a cluster whose removal splits the cluster into two. If we imagine that an electrical 
current is flowing through the cluster, then the "red bonds" would be the first bonds to become red and 
melt, by analogy with a fuse. For loop gases on the cylinder, the previous definition needs to be adapted. 
Let us consider two clusters, percolating from top to bottom and delimited by two defects acting as Peierls 
contours, as in figure |3] (The cylindrical geometry is obtained by identifying the left side with the right 
one.) Flipping one of the red colored w-type faces in figure|3]will create a unique "cluster" encircling the 
cylinder. More generally, a red face is any w-face which, by flipping it to its mirror-state, creates or destroys 
an encircling cluster. In other words, any w-type face formed by quarter-circles coming from both defects, 
is a red face. We will continue to refer to these w-faces as "red bonds" instead of "red faces". 
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The conjectured fractal dimensions discussed in the following are taken from Ill3l . Some of these ex- 
pressions had been proposed before, in Il24ll and others, for the FK representation of Potts models. The hull 
fractal dimension is 



K 



2 — 2Ap'^p±i, dilute phase, 



dh = l + ^=2-2Ao,±i (16) 

'2 — 2Ap,p'±i, dense phase, 

where A^^s is defined as in | [T0) even for r or s non-positive. Beffara ||4| has shown rigorously that SLE^ 
paths have this dimension. The fact that the continuum scaling limit of the hull is actually an SLEk path is 
known rigorously only for a handful of k values. 

The dimension of the external perimeter is conjectured to be 

1 j5 |2 — 2Ap'p±i, dilute phase, 

dep = 1 + ^e(K-i) + -e(i -R) = <^ ' (17) 

|2 — 2Ap±i^p', dense phase, 

where 9 is the Heaviside step function for which 9(0) — j.ln the dilute phase j < k <^, dep = dh = 1 + f , 
which can be understood in terms of CLE^. Indeed, the dilute interval is 2 < k < 4 and, for these values 
of K, it is known that the contours (loops) are almost surely simple, implying that there is almost surely no 
intersection of the contour with itself. This was shown rigorously by Rohde and Schramm 1221 (see also 
Kager and Nienhuis ||T61 ). 

Finally, the fractal dimension of the red bonds is conjectured to be 

i> a (2 — 2Ap'p±2> dilute phase, 

drb = 1+y-^=2-2Ao,±2 = <^ ' (18) 

^ |2 — 2Ap_p'±2) dense phase. 

This fractal dimension is negative on the dilute interval. At first, one might think that (Ts) is amiss, or that 
it does not hold for all ^ < i< < 1 ■ But as explained in Mandelbrot IfTHl , negative fractal dimensions may 
occur in certain physical problems. We will discuss further the issue with negative fractal dimensions in 
section |43] 

We gather in table [T] these predicted values for a selection of dilute models. One of the main goals of 
this paper is to verify some of these predictions. 

3 Measurements of fractal dimensions of dilute models 
3.1 The Minkowski fractal dimension of the defect 

Let S be a subset of II and let e be the mesh of a h5^ercubic lattice drawn on H^. The Minkowski or 

box-coimting definition of the fractal dimension of S is given by 

lnN(e 
e"o Inl/e ' 

where N(e) is the number of boxes that intersect S. In the present case, d — 2 and is replaced by a 
bounded subset of area A. It is natural to introduce a function R(H, V) that measures the linear size of A 
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s = lim (19) 



model (dilute phase) 


K 


P 


c 


dh = dep 


drb 


name, VLM (p , p J 


p 

P' 


—2 cos ? 

K 




1 + f 




percolation, VCMil, 3) 


2 
3 








4 
3 


1 1 

12 


Ising, VCM[3,A] 


3 
4 


1 


1 

2 


1 1 

8 


5 
8 


tricritical Ising, VCM[A, 5) 


4 

5 




7 

10 


7 
5 


19 

40 


3-Potts, VCM[5,G] 


5 
6 


2(\/5 + l) 


4 
5 


1 7 
T2 


^60 


4-Potts, VLM[^,^] 


1 


2 


1 


2 





VCM[5J] 


5 
7 


l(x/5-l) 


1 1 

35 


19 
14 


26 
35 


VCM[?>,5] 


3 
5 


-1 


3 
5 


13 
10 


6 
5 


polymers, VCM[^,T\ 


1 

2 


-2 


-2 


5 
4 


7 
4 



Table 1: Parameters and fractal dimensions for different dilute models. The models with |3 < were not 
simulated in this work but added here for comparison. 

in terms of the mesh e. Here H and V are the numbers of horizontal and vertical boxes in the lattice. The 
simplest definition of R is given implicitly by 

A=(R(H,V)e)^ (20) 

The Minkowski fractal dimension is then ds — lim dl^^^ where 

H,V^oo ^ 

d^^^ = In N (H, V)/ (In R(H, V) - ^ In A) . (21) 

Of course, the numbers H and V, the area A and the mesh e are related. For simplicity, we will use A = 1 . 

The definition pO) of the linear size R(H, V) is natural for studying the defect of £Al(p,p') (see ap- 
pendix[C) as in Il23l . Indeed, (i) each state of a box contains precisely two quarter-circles and (ii) each box of 
the lattice of area A is accessible to the defect. Neither of these two observations holds for the dilute models 
I'£A^(p,p'). Beside the states containing two quarter-circles (w-faces), there are some that contain a line 
segment of length e (v-faces) or only one quarter-circle (u-faces) or nothing at all (the empty face). We shall 
use boundary conditions on the subset S that forbid loops to reach the boundary. In this case, not all states 
are available for the boundary boxes. There are therefore two problems to resolve. 

The first problem is to decide the number of boxes of side e needed to cover each box state. Should 
the u- and v-faces occupy the same area? What about the w-faces? The second problem is to make the 
definition of the area A precise. One could simply decide that each box of the lattice may contain at most 
two quarter-circles so that A — IHVe^, where e is the side length of a box small enough to cover one 
quarter-circle. But one might want to reduce this number of boxes due to the fact that w-faces are not 
allowed along the boundary. Or one can decide to define the area as that corresponding to the number 
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Nmax(e) of e-boxes necessary to cover the largest subset of the lattice that can be occupied by the defect. 
The length of the longest defect, that can be drawn on the lattice, times would be the total area A that 
any observable could possibly occupy, leading to a maximal fractal dimension of 2. Note that the definition 
of the area using Nmax depends on the answer given to the first question. In what follows, we choose this 
interpretation and explore various definitions of R in the case of the hull. 

It is hoped that all reasonable choices will lead to the same fractal dimension for the defects in the limit 
H, V — > oo. Still, some choices might be more appropriate for finite lattices and the rest of this section is 
devoted to see the impact of various choices on the quality of the extrapolation. Note that the following 
discussion is for a strip with the boundary conditions just described. However, it can easily be adapted to 
lattices with cylindrical geometry. 




Figure 4: Covering of the defect for R = ^ ' ^ of a sample 8x8 lattice. 



As a first example of the various weights given to boxes crossed by the defect, consider boxes with only 
one quarter-circle or a straight line segment as half -filled, and boxes with two quarter-circles as completely 
filled. More precisely, if is the area of a tile, then Su = Sv = and Sw = e^, where Si is the area 



occupied by an i-face. A graphical interpretation of this choice is given in figure 4(a) If the numbers 
H and V are even, the defect with maximal Nmax(e) covers the lattice as follows: all boxes not touching 
the boundary are w-faces, while all boundary boxes but one are u- or v-faces. The last boundary box is 
necessarily the empty face. We conclude that, for this choice, the linear size R is given by 

R = (HV-H-V + 3/2)^/^. (22) 

More generally, one may define the two-parameter family 

R^"'^' = (2a(HV-H-V+^)+b)^^^ (23) 
13 




valid when H and V are even, such that = cle , Sv = be and = 2ae , with < a < b < 2a. Then R 



in 1 22 1 is simply ^ ' ^ . Figure 4(b) shows the "area" covered by a defect using this correspondence. Note 
that, unless a = j, e'^ is not in general the area of a box. 

One might also be interested in cases where b > 2q > 0. Now the v-faces cover a larger area than the 
w-faces, and defects with N^ax are found in configurations different from those leading to Ri . We thus 
define 

R^'''"' = (2QV + b(HV-H-2V + 3))'^^. (24) 

This family of coverings contains the choice Su — — Sw = , or R2* ' , depicted in figurels] We 

found this covering to be another relevant choice since it can be argued that a quarter-circle is shorter than 
a straight line segment, and thus that Su should be smaller than Sv Because the range of validity of | [23) 
and | [24| do not overlap, we might as well define 

0<Q<b<2a 

R(a,b)= , (25) 
[r^"-''', b>2a>0 

which is valid as above only when H and V are even numbers, and for the strip only. Although we limited 
our experiments in what follows to members of | [25) , one could in principle extend this definition to include 
the less intuitive situations where < b < a or Sw 7^ 2su. 

Figure |6] shows fits of the same data obtained from various definitions of R. To do these, the numbers 
of u, V and w-faces were stored separately for every configuration measured. We see that all the linear fits 
of dh extrapolate to about 1 .362 at 1 / In R = 0. The exact value is of no importance right now since these 
measurements are imprecise and the lattices used are relatively small. It seems clear, however, that b is the 

^HxV 

important parameter and that, when b = 1, the dh values are closer to their asymptotic values. For 
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Figure 6: Linear fits of the hull dimension dh. against 1 / In R for different covering definitions using the 
model 2?£A^(3,4) on the square geometry. The pair (a, b) refers to R(a,b)- 



example, the two curves with R(i i ] arid R(i i] are far apart, but all those with b — 1 (or b = 1/2) are 
bunched together: this probably means that the v-faces play a significant role in the length of the defect. 
Because the results are closer to their asymptotic values for b = 1 , we choose to work with the covering 

(Q,b) = (l,1). 

3.2 Technical issues 

3.2.1 Simulations on the cylinder 

As explained before, the convergence to the asymptotic value of the fractal dimension of the defect is very 
slow. Despite the improvements made to the algorithm (cf. section [a!2| , lattices of linear size larger than 
H = 512 are difficult to reach. Let us explore the cylindrical geometry using the model 2?£A^(3,4), where 
|3 = 1 . Since no loop counting is required in this case, the algorithm is simpler and thus performs more 
Monte Carlo cycles per second. We measured the fractal dimension of the defect on cylinders with different 
aspect ratios V/H e {1 , 2, 4, 8, 1 6}, where V is the number of boxes along the symmetry axis of the cylinder. 
The entry and exit points of the defect were set in columns corresponding to the azimuthal angle being and 
7t. For each of these ratios V/H, the hull fractal dimensions were obtained for H e {4, 8, 1 6, 32, 64, 128, 256}. 
The plot of linear fits with respect to 1 / In R is shown in figure |7| Obviously, these fits are coarse. Never- 
theless, it is already clear that they converge to a narrow window around the predicted value of equation 
| (T6) , dh = = 1 .375 (horizontal line), as does the curve corresponding to the measurements on the square 
strip. Still, none seems to aim right at 1 .375. Furthermore, the cylindrical geometry might not speed up the 
convergence towards the asymptotic value dh- 

These results call for some discussion. First, the results for the two geometries, that of the strip and of the 
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Figure 7: Linear fits of dh as a function of 1 / In R for cylinders of different V/H ratio and for the square. 
The horizontal line marks dh = 11 /8 for VCM (3,4). 



. — -HxV 

cylinder, behave quite differently. On the longest cylinder V/H = 1 6, the dh 's are significantly higher 
than dh = 1 .375, as opposed to those of the square geometry which are all lower, except for H = 4. Second, 
boundary effects cannot be the only reason for the difference between the two geometries. As discussed 
3.2.2[ significant changes in dh are seen close to the boundary only in a region whose width 



in section 



is a fraction of H. Therefore, the long cylinders V/H =4,8 and 1 6 should have approximately the same 

^HxV 

dh 's. What we see is probably the effect of configurations with larger and larger winding number. 
The defect can indeed wrap aroimd the cylinder as it progresses along it. Intuitively, one would think 
that defects winding, say, twice around the cylinder are longer in average than those with no winding. 
Because the number of configurations with large winding numbers increases with the ratio V/H, a longer 

— ^HxV 

cylinder will have defects with larger dh . Third, one might be tempted to do the measurements of dh 
on cylinders with V/H — 2, because the measurements for coarse meshes are slightly closer to the predicted 
value and the slope for larger meshes is small. Of course, if the choice for these cylinders was based on this 
reason, this would be cheating and it is not clear that the fit over the whole range of R would be any better 
since it would also require non-linear terms. 

3.2.2 Boundary effects and the advantage of the cylindrical geometry 

We now explore the effects of the boundary on the measurements of observables. As the lines below will 
show, it is possible that, in a certain limit, dh might actually be a function of the distance from the boundary. 
We shall try to give a proper definition of the fractal dimension as a function of this distance, but let us first 
consider the effect of joining two lattices with different fractal dimensions. Suppose that two square lattices 
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of the same linear size H are covered by different models, VCM (p , p ' ) and VCM ( q , q ' ] such that ^ 
whose defects have dimension dm and dhi, respectively, with dhi > dhi- And consider configurations on 
a new lattice, obtained by joining the two above. Which (3 should give the weight of the loops that overlap 
the two sublattices? This is a subtle problem and, for the sake of the present argument, it might be enough 
to confine loops to either sublattices. For sufficiently large V and H, the part of the defect in the first half 
will be of length ci H'^'^' and ciH'^'^^ in its second half, for some constants Ci and Cz- Thus, its total length 
is Ci H'''^' (l + H''^'^-^^'^'^' 'ci/ci ) and the fractal dimension in the whole region is 

ln[ciH'i'^'(l+H''^'^^-'i'^i'c2/c,)l/lnH — > dhi, (26) 

where R is approximated by H. Note that this argument can easily be modified to describe two lattices of 
unequal sizes. This shows a simple fact about fractal dimensions: if a geometric object S is studied on a 
region D c R'^ of a finite lattice, then its fractal dimension dg over D will be equal to dg"™, that is to the 
fractal dimension over the subregion Emax C D where it is maximal. 

The present situation is slightly more complicated as dh might be varying continuously with the dis- 
tance from the boundary. What does this mean? The definition |T9l l applies to a subset S C or to a 
subset S of the cylinder as a whole. To define a "local" fractal dimension, consider a cylinder of fixed ratio 
r — length/perimeter. Let L be the length of the cylinder and < l] < I2 < L and let C12 be the annulus 
along (or band around) the cylinder containing points at a distance I from one extremity with li < I < I2. 
Then one can measure dh on C12 by replacing R'^ with C12 in the discussion leading to (19) . Numerical 
estimates of dh can be obtained by covering the cylinder with larger and larger H x V lattices, with r — V/H, 
and counting the number of intersections of the defect in a given Ci 2. Note that the defect can meander out 
of a given annulus, come back to it, leave again, and so on; all of its intersections with C12 must be counted. 
Unfortunately this fractal dimension dh — dh( 1-1,1-2) is difficult to measure numerically. The narrowest 
annulus on the small cylinder H x V = 8 x 8 already takes up one eighth of the cylinder's length. Since the 
fractal dimension is obtained as an extrapolation over several lattice sizes, the smallest section of the cylin- 
der where this dh(li , I-2) can be measured effectively is of length I2 — l-i — 1-/8 — (Ve)/8. Unfortunately, 
some quick explorations have shown to us that a circle at a distance L/8 from the boundary is almost free 
of the boundary effects. The following experiment allows to probe boundary effects. 

Let C be a circle at a fixed distance from the boundary of the cylinder. This circle C has (standard 
or fractal) dimension one. We consider the set c C C of intersections of the defect with C. The fractal 
dimension d^ of this set is therefore < dc < dc = 1 . It is this definition that is used in figure |8] where dc is 
measured as a function of the distance from the nearest extremity, in units of He. Cylinders with (H, V) = 
(32,512), (256,512), (64, 128) and (64, 1024) were used. All the rows of the cylinders were measured, but 
only those at a distance < I < He are shown. Moreover, only a subset of the available data is drawn; the 
size of the datasets for 64 x 1 28 and 64 x 1 024 was reduced by a factor of 2 while the one for 256 x 51 2 was 
by 8. Two observations clearly stand out from figure |8] the fractal dimension seems to "feel" the boundary 
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Figure 8: Fractal dimension of the intersections between the defect and circles C on cylinders 
(H,V) = (32,512), (256,512), (64, 128) and (64,1024) for VCM[3,A). 



roughly up to a distance of l/(He) = 1 and this distance is independent of both the ratio V/H of the cylinder 
and its mesh size e. The following simple argument relates dc to the fractal dimension dh of the hull. 

Let N (H, V) be the number of intersections of the defect with H rows chosen far away from the bound- 



aries of a long cylinder. (We use the notation of section 3.1 ) Choose the circle C to be one of these H rows. 
On average, there will be N/H intersections with this circle C and the fractal dimension dc introduced 
above is the limit of 

ln(N(H,V)/H) 
InRc(H) 

as H ^ oo. As before, Rc(H) stands for the number of "boxes" necessary to cover C. We choose this 
number to be the maximal number of intersections with C, namely 2H — 1 . However, dh is the limit of 
In N (H, V)/ In Rh (H, V) where now Rh (H, V) is the square root of the maximal number of intersections of 
the defect with the H rows, that is approximately VJhP. So, for large H's, dh ~ In N / In H and 

_ ln(N(H,V)/H) lnN_ _ k 

ln(2H-l) ^InH 2' 

Since dh G [1 , 2], the dimension dc lies in [0, 1] as desired and, for VCM. (3, 4), dc should be |, very close to 
the value 0.37 seen in figure|8] 

This result suggests yet another experiment; one in which dh is measured on the annulus Cbuik which is 
defined for li /L = H/(3V) and I2/L = 1 — li /L, that is, an annulus excluding the region where boundary 
effects are felt. Figure |9] shows the data for the measurements of dh on Cbuik for cylinders with ratios 
V/H G {1,2, 4, 8}. There are four other sets of points all converging to the same value, around 1 .35. These 
are the measured dh on the complement Cboundary of Cbuik in the cylinder. Results remain preliminary as 
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they are limited to H < 256 only. But they are striking. The finite fractal dimensions on Cbuik all converge 
quite close to the expected ^ and all those on the complement to another lower value. These results do 
not allow us to determine whether the fractal dimension on Cboundary and that on Cbuik are distinct. Indeed, 
in figure |8j the points for the largest lattice (256 x 512) clearly stand over those of the other lattices when 
the circle C is chosen close to the boundary. They seem to indicate that the limit of dc is not yet reached. 
But even if these two fractal dimensions on Cboundary and Cbuik were the same, these results do allow to 
conclude that the rate of convergence to their asymptotic value is different. This is the main conclusion of 
this long analysis. 
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Figure 9: Measurements of Cboundary and Cbuik on cylinders with V/H — 1 , 2, 4 and 8. 



Figure |9] shows how to measure efficiently any ds in the bulk whenever the rate of convergence of 
d^^^(l) depends on the distance I from the boundary. If we assume that the function d^^^(l) is ap- 
proximately conformally invariant, then the effects of the boundary will be definitely more important for 
measurements on the square. To see this, map conformally the cylinder onto the disk by the exponential. 
Because the circumference of the cylinder has H sites, the annulus that is at distance H/3 from the extrem- 
ity will be sent onto a circle with radius r ~ e^^"/^ Tdisk < 0.125 Tdisk- The Schwarz-Christoffel map that 
sends the disk onto the square will change slightly the form of this inner circle close to the center, but it 
will not change the fact that a minute number of boxes of the square geometry actually lie in the bulk (see 
Langlands, Lewis, and Saint- Aubin flTl ), approximately one hundredth of the total number of boxes. This 
precludes obtaining a good measurement of ds in the bulk using square geometries. Indeed, even for a 
256 X 256 lattice, one would have to count the intersections of the defect with a small square at the center 
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of the lattice of about 30 x 30 boxes. 

We therefore propose to restrict our study to the fractal dimension ds in the bulk, and use the cylindrical 
geometry to measure it, with cylinders of ratio V/H — 2. For simplicity and to minimize boundary effects 
in our experiments, we define the bulk as the region enclosed by the annuli at distance li — W/2 and 



I2 — 3H/2 from one of the extremities of the cylinder To extend the discussion in section 3.1 we must seek 
the maximum number Nmax of intersections the defect can have with this bulk section. Again we assume 
that faces u and v count for one intersection and w for two. For (Q,b) — (1,1), the most dense configuration 
has H(2H — 1 ) such intersections and we define R to be the square root of this number. This choice of the 
size ratio V/H — 2 and bulk region allows to use half of the boxes for measurements with only a reasonable 
increase in the total number of boxes. This seems to be a good compromise. 

3.2.3 The distribution of the winding number 




<p = -n 




a = l — 









-71 TT 

(a) Winding number a = 0. (b) Winding number a = 1 . 

Figure 10: Sample configurations with winding numbers a = and oc = 1 . 



Recall that the entry and exit points of the defect are chosen to be at locations corresponding to azimuthal 
angles equals to and n, respectively. We define the winding number a of the defect by 

7T — (p 



2n 



(27) 



where cp is the azimuthal angle of the defect. A configuration will be said to have winding number oc = 
if the angle of the defect q) has increased by n while going from one extremity to the other. Note that, with 
this definition, a variation of — tt of the angle would amount to a winding number <x — \. Two sample 



configurations having winding numbers and 1 are shown in figure 10 The definition (27 1 is asymmetric, 
and rather unpleasant. However, as will be seen soon, this asymmetry is actually welcome. The simulations 



reported in the previous paragraphs 3.2.1 and 3.2.2 were done starting the thermalization with a defect 
with zero winding number. This has the consequence that no configurations with odd winding numbers 



are considered in the Monte Carlo integration. Indeed, the upgrade algorithm (section[A.l[ always changes 



the winding number by zero or two units. For example, in figure 11 the flip of the gray box changes the 
winding number from to 2. 
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Figure 1 1 : The algorithm preserves the parity of the winding number. 



The exclusion of all configurations with an odd winding number is a serious problem. One could try to 
overcome it by measuring separately the distributions of the defect's length for configurations with even 
and for those with odd winding numbers. For the latter, it would amoimt to start with a configuration in this 
set. But then the determination of the distribution of the defect length for the whole space of configurations 
would raise the question of the relative probability of the two sets and this seems extremely difficult. As 
it turns out, this problem does not occur! Indeed, the probability distribution of the winding number is 
obviously invariant through a mirror containing the cylinder axis, corresponding to cp — cp. Due to our 
definition, however, this symmetry sends a configuration with winding number onto a configuration with 
wuiduig number 1 . More generally, configurations with wuidrng number a are sent onto configurations 
with winding number 1 — a and are therefore equiprobable. Note that, if the entry and exit points had been 
put in the same column, that is with the same azunuthal angle, the problem of discarding the odd or the 
even configurations would have been a serious one. 
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Figure 12: The discrete probability density of the winding number on the V/H — 1 6 cylinder for VCAi (3, 4). 



Figure 12 shows the probability density function of each winding number within its set, either even 
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(square dots) or odd (disks), measured on the cylinder with ratio V/H — 1 6 and H = 8 for the model 
VCM (3, 4). It supports quite nicely the previous observation, as the two histograms are mapped onto one 
another by the symmetry a 1 — a. It is interesting to note also that quite large winding numbers are 
reached. The histogram box for 10 is readable, but larger winding numbers were also reached. Pinson 
||2T] and Arguin Q have considered the probability of Fortuin-Kasteleyn configurations having a given 
homotopy on the torus. Their results may be applicable to long cylinders or, at least, give a reasonable 
prediction for this geometry. 

4 Results 

To estimate the fractal dimension of the hull for a given lattice size, we took m measurements, usually with 
m = 1 000, of the length L of the defect in the bulk region of the V/H — 2 cylinder. Each pair of consecutive 
measurements is separated by a fixed number A of Monte-Carlo iterations which varies according to the 
size of the lattice and the model itself. This procedure was carried out simultaneously on n machines, with 
n > 20, and for cylinders with H = 8, 16,32,64, 128,256,512. The average fractal dimension of a given 
cylinder of size H is then obtained by averaging the defect's length N (H, V) = L of equation < [2T) over all 
the Ta X n data: 

1 ^ ^ 1 T 

E^iLLi^, (28) 

j=l i=l ^ ' 

where R(H) = {2H^ - H) ^ in the bulk for the cylinders with V/H = 2. The same procedure was used for 
the external perimeter. A t5^ical dataset obtained by this method is shown in table |2] 



H 


8 


16 


32 


64 


128 


256 


512 


dn 


1.2772 6 


1.2923 4 


1.3045 4 


1.3143 4 


1.3215 5 


1.3279 7 


1.3323 4 



Table 2: Hull fractal dimension measurements for P£A^(3,4) on cylinders H = 8, 16,32,64, 128,256,512 
with V/H — 2. The notation 1.2772|6 means that the 95% confidence interval's half width is 0.0006 and 
centered on 1 .2772. 

The extrapolated fractal dimensions of the hull and external perimeter for all models studied are shown 
in table|3]with their 95% confidence interval. The meaning of our confidence interval is as follows: if the 
experiment was repeated many times, then an average of 5 fits out of 1 00 would predict an average value 
d at 1 / In R — > that lies outside the 95% confidence interval of table [s] We see that the extrapolations 
are in good agreement with the theoretical values given by equations | [T6) and (TT) . Still, some of the 
confidence intervals do not contain the theoretical values; this is especially true for the external perimeter 
in VCAi (1,1). Even though the definition of the confidence interval does not imply that it should contain 
the theoretical value with 95% probability, one would like to better understand these small departures. 
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model name, VCMip^p') 


K 


13 


dh = dep 




dep 


percolation, VCMil, 3) 
Ising, VCM(3,4) 
tricritical Ising, VCM (4, 5) 

3- Potts, VCM[5,6) 

4- Potts, V£M[],]) 

VCM[5J) 


2 
3 

3 
4 

4 
5 

5 
6 

1 

5 
7 



1 

Vi 

l(^ + ^/5) 
2 

i(V5-l) 


f ~ 1.3333 
^ = 1.3750 

I = 1 .4000 

II ~ 1.4167 
f = 1 .5000 
11- 1.3571 


1.3328 13 
1.3757 13 
1.4025 14 
1.4197 11 
1.5021 12 
1.3592 20 


1.3332 19 
1.3751 15 
1.3986 19 
1.4119 12 
1.4615 10 
1.3586 26 



Table 3: Extrapolated fractal dimensions for the hull and external perimeter. Fits of the fo rm (3o + 13 1 / In R 
P2/ln^R+ |33/ln^R were considered and model testing led us to set (Si = (see section 



B.2.3 



These discrepancies for the hull and the external perimeter will be discussed in the next subsections. 
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Figure 13: Measurements of the hull's fractal dimension against 1/lnR, and the corresponding fits on the 
V/H — 2 cylinder. The theoretical values are designated by wedges at 1 / In R ^ 0. The pair (p, p ' ) refers to 
VCM{v,v']. 



4.1 The hull 



The measurements and fits of the hull's fractal dimension are shown in figure 13 for all models. The error 
bars are much smaller than the symbols used to depict the data. According to the statistical test used (see 
appendix B.2i, the best fits are obtained for the model Y(R) — (3o + (3i/lnR + ps/ln^ R, where |3i and 
(33 are real constants, R is the linear size of the lattice, and |3o is a real constant to be interpreted as the 
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fractal dimension. Note that the necessity of a quadratic term in 1 / In^ R is always rejected by the statistical 
analysis. We do not have a straightforward explanation for this. The figure underlines the quality of the 



match between theoretical values ( 16 1 and measured ones and also reveals some of the possible difficulties. 

First and foremost, the data show that the dh's, as functions of 1/lnR, are not linear and have large 
slopes. The smallest lattice used in the fits has H = 8. On the interval of H covered by our measurements, 

— ^H=8 

just a little more than half the distance between dh and the theoretical value is covered. More precisely, 
the ratio of the data span |Y(R(512)) — Y(R(8)) | over the distance [dn — Y(R(8)) |, averaged over all fits, 
is only 0.57. Even with data for larger lattices, say H — 1024,2048 and 4096, this ratio would still be 
under 0.7. These extrapolations are obviously difficult. Second, the fits are extremely sensitive to the non- 
linear terms of Y(R). According to table |3] and figure 13 the fits give an extrapolated fractal dimension 



that slightly overshoots the theoretical one, the only exception being the fit for 2?£A^(2, 3). The fit for 
VCAi (3, 4) overshoots very slightly, by about 0.0007, and the others by about 0.002. The latter are precisely 
those whose graphs are the most curved. The departure from the predicted values is even more striking for 
some values of the external perimeter and we shall propose in the next subsection an experiment to better 
understand the systematic error caused by the extrapolation. 

Of course, better estimates of dh would require data for larger lattices. For our algorithm, the major 
obstacles to reach these lattices are the following. First, the smaller the u, v and w weights (see l|7}), the 
slower the simulations. The reason is that, when the algorithm steps on an empty 3x3 block, the probability 
that the block remains empty ranges from 90% for 'DCM{^, 1) to over 99% for VCMiB,?]. Because the 
functions u, v and w increase with |3, the autocorrelation function of the Monte Carlo chain decreases 
faster for larger |3, as there is less emptiness. Second, the dependency of dh on boundary effects favored 
simulations on the cylinder. But accounting for the cyclical boundary conditions slows the algorithm by 
about 5%. Third, the amount of data on the lattices up to H = 256, for (3 ^0,1, are easily stored in the cache 
of the cpus used but, as H reaches 512, more traffic between the cache and the memory is needed as seen by 
the sudden decrease in the number of Monte Carlo cycles performed per second for these larger H values. 



Finally, the model VCAi (5, 7) is the hardest to measure. As seen in figure 13 H = 256 is the largest lattice 
for it. Its fugacity ^{sj) — 0.6180 is the lowest among the models measured, excluding |3(2,3) =0 that has 
its own algorithm anjrway. The u, v and w are small and the cost of creating a loop starts to show, leading 
to rather empty configurations with the problem noted above. 

4.2 The external perimeter 

The definition we used for the external perimeter is inspired by the biased walker of Grossman and Aharony 
IITSi and its variant for loop models introduced in Saint- Aubin et al. 1231 . Starting at one end of the cylinder, 
the walker follows the left side of the defect all the way down to the other extremity of the cylinder and, in 
the meanwhile, counts a unit of length each time it encounters a u- or v-face. When it encounters a w-face 
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made up of two quarter-circles coming from the defect, it still counts a unit of length, but then it does not 
enter the fjord created at that point, unless the two quarter-circles possess azimuthal angles cpi and cpi 
such that |(pi — (P2I = 27T. (See section 3.2.3 for the definition of cpt.) That is, if we interpret the defect as 



a chordal SLE growing from one end of the cylinder to the other, then its azimuthal angle as it progresses 
along the trajectory takes the value q)i when it first crosses the w-face and (p2 the second time. The only 
possible differences of these two values are jcpi — cpil =0 and In. When the difference is In, we consider 
that this w-face is not the entry point of a fjord because otherwise the walker would return to its starting 
point without scanning the defect in its entirety. The external perimeter for two distinct defects is drawn in 
figure |3] 

As observed in ||23ll , there could be more than one definition for the external perimeter in loop models. 
Thus, to confirm the validity of the one chosen, we have compared the extrapolated value obtained for 
the dense phase of 2?£A4(2,3) with that of £A^(2,3), previously measured in Il23l . Since the two models 



belong to the same universality class, as explained in section 2.2 the two extrapolations should be equal. 
For 2?£A4(2, 3) with the above definition, we measured 1 .3368, while one finds 1 .326 for £A^(2, 3] in 1231 . 
Both measures are fairly close to each other and to the theoretical value of dep — 4/3. We thus conclude 
that there is no reason to doubt our definition. 
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Figure 14: Measurements of the external perimeter's fractal dimension, and the corresponding fitted curves, 
for all the previously cited models on the V/H = 2 cylinder. The pair (p,p') refers to 2?£7M(p,p'). 



The measured dimensions dgp for the external perimeter are shown in table [s] Again, we see that 
the agreement with the prediction is generally good, although the results for I?£A^(5, 6), and particularly 
the one for D£M[], 1 ), are farther away from the theoretical ones than the corresponding results of the 



hull. Nevertheless, as visible in figure 14 the fits are good enough to confirm convincingly the theoretical 
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predictions of equation 1 17 1, except maybe for VCM (1,1). The latter model is in the equivalence class with 
central charge c = 1 to which also belongs the 4-Potts model. It is known that, for many models at c = 1 , 
the scaling of geometric objects include logarithmic terms (see for example H)) and thus extrapolation is 
predictably difficult for VCM (1,1). Still, one would like to better understand the systematic undershooting 
for this geometric exponent and, to a lesser extent, for that of the hull. We propose the following experiment 
with this goal in mind. 

To the Monte Carlo data for VCM (1,1) and VCM [5, 6) that we have, we added the point ( 1 / In R, dep ] = 
(0, dg*^"), i.e. the theoretical prediction. (This idea is not new. For example, Asikainen et al. fS) also use fits 
where the theoretical fractal dimensions appear explicitly to determine their error bars.) And we fitted 
the new datasets with this single additional point. According to the F-test (see appendix B.2.3| , it is now 

2. 3 

appropriate to keep both of the terms 1/ In R and 1/ In R, and the coefficients ofY(R) — (3o + Pi/lnR + 

2. 3 

Pi/ln R+ (33/ln R are now all significant. Visually the two fits, with and without the theoretical value. 



are barely distinguishable in the range of measurements as can be seen in figure 15 but they split quickly 
to the left of the datum corresponding to the largest lattice. We therefore propose the following simple 
interpretation for the systematic departure from the predicted values: in the range H < 512, the cubic term 
1 / Ln^ R is large enough to conceal the quadratic one. Without the additional (predicted) point, the F-test 
says that the h5^othesis that (32 is zero cannot be rejected. (Actually, if one includes pi, the dh is similar, 
but the confidence interval is larger) The fit with only Y(R) = |3o + (3i/lnR + (Ss/ln^ R is then excellent. 
But, with the predicted value added, the hypothesis that P2 is zero must be rejected and the new fit goes 
through all the measured dh with equal precision. We made one further check by adding, instead of 
the theoretical value, a set of random points at 1 / In R = 0.05 and 0.1 whose averages sit on the new fits of 



figure 15 with a variance similar to the spread of the data at H = 512. The fit with these random points had 
the same properties as the one with the theoretical value added. The conclusion seems to be that the sizes of 
lattice we used do not allow to measure properly both pi and P3 that describe finite-size corrections from 
H ~ 8 to H ^ 00. The values 1 / In R = 0.05 and 0.1 correspond to lattices of size H ~ 3 x 1 0^ and 2 x 1 C*. 
They are out of reach with our algorithms! 



4.3 The red bonds 



H 


8 


16 


32 


64 


128 


256 


drb 


0.1315 5 


-0.0243 4 


-0.1262 7 


-0.1991 9 


-0.2510 11 


-0.2941 21 



Table 4: The measured fractal dimension of red bonds for the Ising model VCM{3,4) on cylinders 
H = 8, 1 6, 32, 64, 128, 256 with V/H = 2. 



To measure the fractal dimension of red bonds, we use two defects, instead of just one, so that the do- 
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Figure 15: Fits and measurements for the external perimeter of 2?£A4(5,6) and 2?£A^(1, 1) as well as the 
fits with the additional point d*^° at 1 / Ln R = 0. 



main between them is interpreted as a proper "cluster". The fractal dimension is then computed by finding 
the number of w-faces visited by both defects and multiplying this number by two (see section |Z3j. The 
measurement of this observable is difficult since, for models in the dilute phase, not even all configurations 
have at least a single red bond. Moreover, the average number of red bonds is smaller than 1 . Because of 
this, the fractal dimension cannot be estimated using equation | [28) and we used an alternative definition 
that was also used in Saint- Aubrn et al. ||23l : 

1 ^ InNT 



n ^ — 



InR(H)' 



(29) 



where Ni 



^rid Nij is the number of red bonds of the i-th datum in the j-th Markov chain. 



4.1 



The numbers m, n and R(H) have the same meaning as in section ' 

As the data of table |4] attest, there are less and less red bonds in average for each lattice configuration as 
H, V oo. Because of this, we restricted our measurements to 'DCA4(3,4), as it has the most efficient 
algorithm, and confined the lattices to H e [8,256]. The difficulty is well illustrated by the following 
numbers: the average value —0.2941 ± 0.0021 at H = 256 in table [i] required a total of about 10^^ Monte 
Carlo iterations, excluding the number of iterations needed for an adequate thermalization of the startup 
lattices. By comparison, the result 1 .3323 ± 0.0004 at H = 512 for the hull required about the same number. 
In other words, obtaining a measurement for the hull at H = 256 with a similar confidence interval of 
±0.002 would require only about 5.5 x 10^^ iterations, i.e. approximately 1 /1 80 times the number required 
in the red bonds experiment. 

Although some of the measurements are less precise than for the other observables, we obtained the 

27 



drb 




' ' 1/ln R 
0.35 



Figure 16: Measurements of the red bonds' fractal dimension and its corresponding fit, for the dilute Ising 
model VCM (3, 4) on the V/H = 2 cylinder. 



extrapolated drb = —0.621 ± 0.004, in good agreement with the predicted value drb = —0.625 (eq. |[T8|). 
This result is surprisingly close to the theoretical value, considering the extreme extrapolation that had to 
be done: the value at H = 256 remains at about 0.33 from the theoretical prediction, that is 10 times farther 
than the corresponding value was for dep in VCAd (1,1). However, the fit gives a good result because the 
curvature is gentler, and the discrepancy from linearity is rapidly wearing off as 1/lnR decreases. The 



average values as well as the fit for this observable have been plotted in figure 16 

The negative fractal dimension drb is somewhat of a novelty here, but can be easily interpreted. As the 
lattice size grows, the average number of times the two defects touch one another decreases toward zero, 
and the rate of this decrease is measured by the negative fractal dimension drb • 

5 Concluding remarks 

The agreement between the theoretical values and the measured ones (table |3| is very good. It is similar 
in precision to that obtained in Il23ll for the dense models, or even slightly better. Note that the difficulties 
that had to be overcome for the dilute models were quite different from those encountered for the dense 
ones. The improved quality of the geometric exponents obtained here might stem partially from a more 
systematic and finer statistical analysis (see appendix |BJ. 

The proposed upgrade algorithm for the dilute models, and its variants for the models with (3=0 and 
1, was quick enough to provide very precise measurements for lattice sizes uptoHxV — 512x 1024. 
Still, one feels that, unless willing to wait for the next or the second next generation of cpus, the algorithm 
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proposed here has reached its practical limits. New methods, in the direction of those proposed by Deng 
et al. HI, will be necessary to probe further these models. 
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A Upgrade algorithms 

In this appendix, we first present the upgrade algorithm common to all the models. The subsequent sub- 
sections sketch the ideas relevant for particular models. 

A.1 The basic algorithm 

For the dense loop models £A^(p,p'], reviewed in appendix [C) each box is in either one of two faces, 
corresponding to the two w-faces of 'DCA4[p,p'). A straightforward Metropolis-Hastings upgrade can be 
chosen as simply changing one box at a time and checking whether the usual Monte-Carlo condition is 
verified. (See 1231 .) For the dilute loop models P£A^(p,p'), the edges of the nine possible states are not 
necessarily crossed by loop segments, as is shown in figure |T] Indeed, the empty state contains no loop 
segments, two edges of the u- and v-faces are crossed by a loop segment and all four edges of the w-faces 
are crossed. As the loop segment must remain continuous at box interfaces, we could refine the previous 
algorithm by requiring that the edges that are crossed before the change, and only those, remain crossed 
after the change. But then the only effective change would be the flipping of a w-face to its mirror face, 
thus preventing the algorithm from sampling over all configuration space. Because of this, we are forced to 
consider changing many contiguous boxes in a single upgrade step. 

The upgrade step must therefore change an ra x n block of boxes, with m, n > 2. First an ra x n 
block is chosen at random in the H x V lattice on the cylinder Because this block has to fit entirely in 
the lattice, this amounts to placing randomly the upper left box of the block in the first V — n + 1 rows, 
all possibilities being weighted uniformly. Second, the content of the ra x n block is changed for another 
admissible block. To be admissible, the 2(Ta + n) edges of the block must be crossed by a loop if they were 
in the original block, and be free of crossing if the original edge had none. For instance, if the chosen 3x3 



block corresponds to figure 17(a) then an admissible replacement is shown in figure 17(b) while figure 17(c) 



shows a forbidden replacement block. The new block must be chosen uniformly among admissible ones. 
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Finally, the Boltzmann weight u"u-ynv^n„ pN q£ ^j^g configuration where the ra x n block is replaced by 
the new choice is computed and compared with the original one (see l|6)): the Metropolis-Hastings ratio 
decides whether the replacement is to be accepted or rejected. Note that, even though the change is limited 
to the block, the computation of the weight involves counting the number of closed loops (except when 
|3 = 1 ) and might therefore require exploring the configuration at a large distance of the ra x n block under 
consideration. 




(a) The chosen block. (b) An admissible replace- (c) A forbidden replace- 

ment block. ment block. 

Figure 17: The original 3x3 block in (a) with an admissible replacement in (b) and a forbidden one in (c). 



The requirement of uniformity of the block state among admissible ones raises a subtle problem. One 
might think that it is achieved by simply choosing the face of each of the m x n boxes one after the other, 
respecting at each step the conditions on the perimeter This is not the case as the next example shows. 



Suppose that the ra = n = 3 block to be changed is that of figure 17(a) and that the new box states will 
be chosen from left to right, top to bottom. If the first box to be chosen is the upper left one, three choices 
are admissible: the u-face uz that has a single west-north crossing, and both w-faces wi and wz- (As 
2.1.2 the index of the letters u, v and w refers to the order of figure [ij the bottom box being 



in section 



labeled by 1 .) Each of these three faces will be given probability |. It is easy to check that there are two 
possible replacements for the next box, and only one for the last box of the top row. Similarly, two faces 
are admissible for the leftmost box of the second row. So far, twelve different fillings are possible, each 
occurring with probability jj- The difficulty occurs for the box at the center. Two of the twenty-seven 
possible fillings are shown in figure [18] For the filling on the left, two choices are allowed: the U4 or the 
empty face. But for the filling on the right, three are possible, namely the U2, wi and wz faces. This means 
that, starting from the center box, the blocks obtained from the filling on the left will occur with probability 
2^, and the ones obtained from the right filling will get probability j^: this violates the requirement of 
uniformity. To assure uniformity, one has to determine first, for a given block, the number of allowed 
replacements. We found it more efficient to count beforehand these for all the possible edge configurations 
of the perimeter, or boundary state, for the ra x n block, and actually construct a list of the admissible 
replacements. For m = n = 3, there is a total of 11 3 361 possible replacements for the 2048 boundary states. 
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Figure 18: If we chose the face of each box sequentially, then the left block would be more probable than 
the right one. 



Depending on the boundary state, there can be as little as 1 8 admissible replacements or as many as 690. 

We decided to work with 3x3 blocks. The list of possible replacements for larger blocks would take 
up much more memory. Moreover, given some boundary state on the block being changed, the variations 
in probability brought by tentative replacements would likely be in a wider range and lead to a higher 
rejection rate. Such an argument does not hold for 2x2 blocks, the smallest that allow exploration of 
the whole configuration space, for which the list of replacements is short and the acceptance rate is high. 
However, some experimentation with this shorter list shows that the autocorrelation between upgrades, 
and therefore the number of upgrades between statistically independent measurements, is very high. We 
found the 3x3 blocks to be a good compromise. 

A.2 Improvements 

The Ising model is described by the dilute phase of 'DCAi[3,4). Since the loop fugacity of this model 
is |3 — 1 , there is no need to know the number of loops crossing the m x n block when computing the 
Boltzmann weight. This simplification greatly speeds up the algorithm because loop counting is by far its 
most time-consuming component. 

It is more difficult to improve the algorithm for models with (3 G (0, 2] \ {!}. The rest of this paragraph 



is devoted to this problem, while the case (3=0 will be discussed in subsection A. 3 One can distinguish 
between models according to whether 0<(3<1or1 <|3<2, because models of the first type will tend 
to be filled with less loops than those of the second type. Most configurations of both types of models have 
large loops, some being of length of the same order as that of the defect. 

Our first attempt to compute Boltzmann weights was rather naive. We simply followed every connected 
path that crosses the m x n block until it returns to its starting point or, if the path is part of the defect, 
reaches the boundary of the lattice; that works, but it is slow, especially for large lattices. It is slow due to 
the presence of the defect and, potentially, of long loops. The problem is particularly acute for those models 
with large |3, e.g. the dilute phase of 2?£A^(1 , 1 ), because the predicted hull fractal dimension of the defect 
is large and lots of loops are present. 

A way of reducing the impact of the defect crossing the m x n block is to assign a time order to each 
edge it crosses, starting from its entry point down to its exiting point. Because every edge crossed by the 

31 



defect is then identified in some way as belonging to it, it is possible to distinguish the defect from a loop 
during loop counting. Moreover, the time ordering allows to identify when the defect enters first the block 
and when it leaves it for good. Since these two edges cannot change when going from the original to the 
replacement block, it allows for quick identification of those situations where the defect generates a loop or 
absorbs one intersecting the block. 

Large loops are common, some with size commensurate to that of the defect. When the m x n block 
selected for the upgrade is crossed by one of these, a slowing down of the algorithm similar to that en- 
coimtered for the defect is observed. To get over this problem, we time-ordered the edges of each loop, 
starting from an arbitrary edge on its path, and also assigned a unique number to each loop. During the 
loop counting phase of the former block, the number of different loops crossing the block is now obtained 
very quickly: it is simply the total of different loop numbers. The time order is useful because it allows the 
algorithm to know how the different edges of the same loop or defect are connected to each other outside 
the ra X n block. Like for the entry and exit points of the defect, those external connections will not change 
after replacement of the ra x n block. So, with this information, during the loop coimting step for the re- 
placement block, the algorithm does not have to follow the path of the loops (or defect) which lies outside 
the block since the reentry point is now known. 

A.3 The case (3=0 

The percolation model, as described by the dilute phase of X'£A^(2,3], needs a special algorithm to be 
simulated efficiently since its |3 = loop fugacity forbids the presence of loops. It is still possible to choose 
an m X n replacement block on the sole basis that it suits the boundary conditions of the block, but then if 
the chosen configuration creates a loop, this choice will have to be rejected. To get rid of this problem, we 



need to know what the external connections of the block are, as described in section A. 2 this information 
is easily retrieved if the defect is time-ordered. All that is left to do is to choose, with uniform probability, 
a replacement amongst the blocks respecting both the boundary conditions and external connections. In 
figure [19| an example of a possible replacement admitting no loop is shown for a 3 x 3 block on a lattice 
with cylindrical geometry. 

To obtain a new list of replacement blocks appropriate for this model, the external connections of the 
m X n block have to be taken into account, so it is necessary to be able to enumerate them all. To do so, 
one may use a diagrammatic representation, an example of which can be seen in figure |20] This example 



of a diagram corresponds to the external connections of figures 19(a) and 19(b) Under this representation, 
the connected edges are linked pairwise by arches, and the entry and exit edges of the defect in the block, 
depicted by pale dots in the figures, are linked at infinity by vertical lines. The linking arches may not cross 
each other, but they may go "under" the vertical lines, that is without intersecting them, as the defect may 
wind around the cylinder without self-intersecting. For given block boundary conditions, consisting of p 
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(a) The block to replace. 



(b) An admissible replacement. 



Figure 19: A possible transition for a 3 x 3 block on the cylinder yielding no loop. The first and last boundary 
edges crossed by the defect are marked by pale dots, while the remaining connected edges are marked by 
black ones. 



Figure 20: Diagrammatic representation of the external connections shown in figure 19 
aligned by starting from the top left edge and following the boundary clockwise. 



The points are 



connected and 2m + 2n — p imconnected edges, there are 

Ep = ^P(p-nCp_2, (30) 

different external connections, where Ci, = (1/2) ~ (i/2-i)' obtain the new list of replacements, one 
must first fix the boundary conditions and external connections of the ra x n block and then, using the list 



of section A.l check whether any of the suggested replacements yields loops; those that do not generate 
loops belong to the new list. This testing has to be repeated for all possible external connections matching 
these fixed boundary conditions. Finally, repeating this process for all 2048 boundary conditions completes 



the list of replacements. For a 3 x 3 block, we find, using equation 1 30 1, a total of 

N = ^ (\^)^^ " 144408 



i=2 
i even 



different external connections, leading to a list of over 3.8 x 10^ replacements. This is to be compared with 
the 113 361 possible replacements of the standard |3 > list. 

When the (3=0 algorithm steps on an empty m x n block, nothing can be done since loops are not 
allowed, and so the algorithm skips this block and chooses a new one. According to equation 0, the |3 = 0, 
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or A = 7t/8, model is the most diluted of all the |3 > models. That is, the ratios of the u, v and w weights 
with that of the empty face are at their smallest. Empty blocks therefore occur frequently for (3=0 and 
the present algorithm takes advantage of this. Even though the |3 = algorithm is involved and its list of 
replacements is heavy, it is the second fastest algorithm, second only to the one for |3 = 1 . By second fastest 
we mean that the algorithm for (3 = computes in average the second most Metropolis-Hastings iterations 
per second. 

B Statistics 

In this section, we explain the procedures we used to obtain warm-up intervals, the statistical formulas we 
used to obtain the confidence intervals on measurements, and the extrapolation procedure. 

B.l Warm-up interval 

The warm-up interval, or burn-in period, is the average number of Monte Carlo cycles needed to attain ther- 
malization. To find a reliable warm-up interval for each of the models I?£7W(p,p'), one may use the stan- 
dard procedure (see Fishman 1141 for instance). From any H x V lattice configuration, start the Monte Carlo 
algorithm for n independent Markov chains and take m measures of the length L of the defect in the bulk. 
Each measure is to be separated from the next one by A Metropolis-Hastings (MH) cycles. Finally, plot 
the measurements, averaged on these n chains, against m: the abscissa where the average values stabilize 
provides a warm-up interval. 

The approximate point of thermalization for VCAi (3, 4) on the cylinder of size H = 51 2 and V — 1 024 is 
represented in figure |2T|by a vertical line at m = 400. After that point, the average value fluctuates around 
L = 6500 which, according to equation 
This is to be compared with the final measurement of table |2j dh = 1 .3323 ± 0.0004. To illustrate 
the reliability of this thermalization point, we plotted two datasets in figure [21] the dark one corresponds 
to a set of measures started from a configuration having the smallest possible defect length in the bulk, 
that is L = Lmin = 512, while the light one corresponds to a set started from the maximal defect length 
L = Uax = R^ =523 776. 

B.2 Statistical analysis 

We give the details here of the methods used for obtaining confidence intervals on the different measures, 
on how the linear regressions, or fits, of the data were done, and also what processes were used to discrim- 
inate between a good and a bad linear regression. For more details on these subjects, see Draper and Smith 

a. 



211 with R = V523 776 ~ 723.7, corresponds to dj^''^ = 1.33. 



34 




2000 - 



200 400 600 800 1000 

Figure 21: Average defect length L in the bulk as a function of the number of measures for a VCM (3, 4) 
512 X 1024 lattice on the V/H = 2 cylinder. This plot was realized for n = 20, m = 1000 and A = 200 x 10''. 



B.2.1 Confidence interval 

For an experiment consisting of n independent Markov chains, e.g. n computers or processes, and a total 
of m measurements Qt j per chain, where i labels the chain and j the datum, the unbiased estimator of the 
expected value E[Q] = Q of an observable Q is given by the grand sample average 

i=l 

where Qt = -j^ ^j=i Qi^-i average value of the i-th Markov chain. To obtain a confidence interval on 

Q, we first need to compute the unbiased sample variance 

1 



n- 1 

1=1 



In this work, some of the observables d^^^ were measured with a small number n of Markov chains 
(n ~ 20). In these cases, it is better to replace the approximate 95% confidence interval 2(t/ -/ti by the correct 

(31) 



Q-Tn-i(0.95)-^, Q+Tn-i(0.95)-^ 



where Tn-i is the inverse Student t-distribution with n — 1 degrees of freedom. For n = 20 observations, 
the factor T19 of equation | [3T| | equals approximately 2.093. 
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B.2.2 Linear regression 

The estimated fractal dimension of the different observables in the continuum scaling limit, as R ^ oo in 
| [2T| , is obtained by extrapolating the linear regression of the measurement data acquired at finite R's. For 
an experiment of ni + n2 + • • • + = n total measures (sample points) of the observable, with the 
number of measures for the i-th lattice size (e.g. Hi — 32, Ht+i = 64, etc.), the most general linear model 
with p parameters is: 

Yi = po + (3iXii +|32Xi2 + --- + |3(p_i,X|(p_i,+ei, (32) 

where Y|, 1 = 1 , 2, . . . , n, are the sample points, pk, k = 0, 1 , . . . , p — 1 , are parameters associated with the 
p independent variables X^k, and ei ~ N (0, cr?) is the i-th random error associated with Y^. We shall be 
interested in pol3momial fits, for which Xij = X|, with Xt = 1 / In R(Hi, Vi). Moreover, Y^ corresponds to the 
measured fractal dimension d^' at R(Hi, Vi). The equation | [32) can be rewritten in matrix notation as 

Y = Xp + e, (33) 

where Y and e are n x 1 vectors, |3 is a p x 1 vector, and X is a n x p matrix. Note that the first colvunn of 
X, related to (3o, is filled with Ts. The fitted equation, 

Y = Xp, (34) 

is such that Y = E[Y] and p = E[p] are unbiased estimators of Y and p, respectively; these are the quantities 
we need to evaluate. 

The difficulty with the present datasets is that the observables d^' have not been measured with the 
same precision, that is, their variances depend on the size of the lattice. And when Var[e] = cr^V ^ cr^ln/ 
with cr^Vij — Cov[ei, ej] and a a positive constant, one cannot rely on ordinary least squares (OLS) for 
obtaining an unbiased linear regression, but must count instead on weighted least squares (WLS). In other 
words, if the sample variances o^Vu = Var[ei] — Var[Yi] in | [32) cannot be considered constant throughout 
the dataset, one cannot use OLS. Moreover, the WLS method may be used only if V is diagonal, i.e. if the 
different are uncorrelated. Fortunately, this is the case here as our Markov chains are independent. 

The idea behind WLS is to multiply both sides of equation | [33) by the constant n x n matrix P — V^^ 
in such a way that the variance of the transformed error Pe becomes constant amongst the dataset; indeed, 
Var[Pe] — cr^ln since P is symmetric. In practice, one may set cr^ = 1 since this constant is the desired 
variance of the transformed variable PY, in which case V — diagjaT^, ctJ^, . . . , ov^^}. The expressions for 
P, Var [p] and Var [Y] in the WLS method are obtained by replacing X with PX and Y with PY in the usual 
OLS expressions. The WLS expressions are then 

P = (XW-'X)"^ XW-^Y 
Var[p] =(j^(XTV-iX)"' (35) 

Var[Y] = cT2xVar[p]XT. 
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Since we are interested in polynomial models, the last expression of 1 35 1 may be rewritten to yield Var [Y] 
as a function of x: 



Var[Yx] = cr2xTVar[p]x 

= CT^ (1 , X, x^ . . . , xP-' ) Var [p] (1 , x, x^, . . . , x^-^ )^ ^ E(x), (36) 

where = (3 is the linear regression. Using | [36) , the 95% confidence limits of as a function of x is 



Y,-Tn-p(0.95)E(x), Y,+T^_p(0.95)E(x) , (37) 



with Tn-p (0.95) as in ||3T). 




l/lnR 



Figure 22: The sample points, the linear regression and its 95% confidence limits for dh of I?£A^(3,4] on 
the V/H = 2 cylinder. 



In figure 22 the results of fitting the model Yt = (3o + Pi Xt + PsXi to the dataset of the observable dh 
(eq. ^) for V£M(3,4) on the V/H = 2 cylinder are shown. The average of dh for each run at each size 
H X V is shown by a small dot. The bottom and top curves are the confidence limits obtained by | [37) , and 
the central curve corresponds to the fit Y^. The 95% confidence interval at 1 /InR = is [1 .3744, 1 .3770], thus 

= 1 .3757 ± 0.001 3 (cf table|3). 

B.2.3 Model testing 

Once a linear regression has been obtained, it is necessary to attest its validity and quality. To do this, the 
popular Pearson correlation coefficient ^[Li (Y^i ^ Y)^^ ^r=i (Yt — Y)^ can be of help, but alone may lead 
to biased results, as discussed in § 11.2 of Draper and Smith ||9j, especially if an extrapolation from the fit is 
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intended. To complement the correlation coefficient, one may use the ¥-test that compares different models 
of linear regression for Y. To do so, we must first fix the full model that contains all the |3 parameters that 
are sensible to use. For instance, our linear regressions used a polynomial model of order 3. Higher order 
polynomials would have required data on more lattice sizes. The full model is thus 

Yi = |3o + Pi Xi + PaX? + 133X3 + a. (38) 

Different models are then compared with the full model. The F-test is used to verify linear hypotheses. 
For example, we may want to compare the quality of the fit obtained, say for a cubic model like | [38) but 
with pi = 0, or an even model where pi = Pa = 0. We might also want to test a hypothesis of the form 
po + 2Pi =4 and po + Pi + Pa = — 1/ where there are now two relations to be satisfied at the same time. 
More generally, a linear hypothesis can be written as 

Ho : R(3 = r, (39) 

where R is an ra x p matrix providing ra linear relations amongst the p's, where q of these restrictions are 
linearly independent. The m x 1 vector r contains the constants of the m relations. 

To verify such a hypothesis, both the estimator pr of the restricted model, Y = XPr + £/ 



-1 



R(XTV-^X) RT (r-R'p), (40) 

and the estimator (3 of the full model are computed. Then the residual sum of squares SSE(pr) and SSE((3) 
for both models are obtained. Their estimators are defined by 

SSE(p) = (Y-Xp)^(Y-Xp), (41) 

and similarly for p^. Finally, the F-test for the hypothesis | [39) consists in comparing the ratio 

/ SSE(p;)-SSE(p) \ /n-p\ 

SSE(p^ jl^J' ^''^ 

with the value z for which F(q, n — p, x) dx — A, where A e [0, 1] (we chose A = 0.95 in this work) and the 
^-distribution is 

1 \"^i/2 / ^ X -(Tt,+n2)/2 

nr^^^^2^-)-W7^i-^ x^-ifl + ^x") , (43) 



with B(x, y ) the Beta function. Now, if < z then, with probability A, the dataset does not provide suf- 
ficient proof that Hq has to be rejected. On the other hand, if > z, then the dataset shows that the 
hypothesis Ho is implausible and should be rejected. Different scores coming from different hypotheses 
concerning the same full model may be compared with each other, the lowest score being the most plau- 
sible. We use this test to check whether one or more of the coefficients pi could be set to zero. When this 
is a possibility, we note that the estimated fractal dimension obtained from the restricted model is, most of 
the time, closer to the theoretical value than that of the full one, and its confidence interval is also smaller 
For the dataset of the d^""^ ot VCM (3, 4), the full cubic model of l|38) gives = 1 .377 ± 0.007, while it is 
1 .3757 ± 0.0013 for the restricted model having the lowest score, that is the one with p2 — 0. 
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C Logarithmic minimal models 



0™ 



Figure 23: The two equally weighted faces of CM. 



Here the lattice models used by Pearce et al. 11201 to define the logarithmic minimal models £Al(p,p') 
are recalled. They were used in 1231 for simulations of the dense models. The underlying loop gas is based 



on two elementary faces, illustrated in figure 23 At the isotropic point, the two faces are given an equal 



weight pi = P2 which can be fixed to 1 . The partition function of the loop gas is then 



(44) 



where C is the set of all possible loop configurations, |3 the loop fugacity, and N the number of loops in 
a given configuration. These models do not correspond to a particular solution of the VCM. models (l|6) 
and Q) as no value of A e [— 7t/2, 7t/2] in any of the critical regimes found in Blote and Nienhuis 151 yields 
u = v = 0, w^O and no empty face. The simulations in ll23l were done on a cylinder. The boundary 



conditions consisted in half-circles added at the extremities of the cylinder, as shown in figure 24 Note 
that, with these boimdary conditions, there are \C\ =2'^^^ possible configurations. 
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Figure 24: An example of an 8 x 8 configuration on the cylinder. 



It is these lattice models whose continuum scaling limit was given the name of logarithmic minimal 
models £A^(p,p'). Their loop fugacity is given by (3 — — 2cos (?) with R = p'/p and they are believed to 
be described by logarithmic CFTs whose central charge and conformal weights are given by l|9| and (TO) . 
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